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Abstract 



The Kirchhoff elastic theory of thin filaments with spontaneous curvature is 
employed in the understanding of the onset of the kink transitions observed 
in short DNA rings. Dynamical analysis shows that when its actual curvature 
is less than some threshold value determined by the spontaneous curvature, 
a circular DNA will begin to buckle to other shapes. The observable and the 
dominant deformation modes are also determined by dynamical instability 
analysis, and the different effects of Zn 2+ and Mg 2+ ions on DNA configura- 
tional properties are qualitatively discussed. 
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This paper is motivated mainly by a recent experiment on short DNA rings.!'! In this 
experiment, Han et al. synthesized a kind of DNA sequences known to have intrinsic bending 
tendency (in its undistorted state this kind of sequences will form an axial bend of at 
least 30° per helical turn.!) They used these sequences to investigate whether DNA can 
actually be kinked (forming large bendings over only few base pairs) and, if kinks really 
turn up, what are the necessary conditions to produce them.!H Kink deformation has been 
theoretically proposed to be a very important mechanism for wrapping DNA around the 
nucleosome particles ever since the seventies,! - ! and a recent x-ray crystalline structure 
of nucleosome core particle also demonstrated that DNA is not uniformly bent but has 
maximal and minimal curvatures at different positions.! However, there had been no direct 
evidence for the existence of kinks in DNA before the experiment of Han et a/.!'! which 
clearly demonstrate this. To their great surprise, they found that kink formation is closely 
related to total chain length and solvent ionic conditions. For DNA circles with 168 base 
pairs (bps) or 16 helical turns, the chain is often kinked and form polygonal shapes if the 
solution has an appropriate Zn 2+ or Zn 2+ /Mg 2+ concentrations. However, kink transition 
does not take place in Mg 2+ solution alone or in Zn 2+ /Mg 2+ solution with high concentration 
of Mg 2+ ions. For DNA circles of 126-bp or 12 helical turns, no kink will ever occur, no 
matter what kind of solvent ionic conditions.!'! Trying to interpret their novel observations, 
Han et al. suggested that axial stresses exist in 168-bp DNAs and they are the main reason 
for the observed kink deformations. 

In this article we attempt to check, from a theoretical point of view, the validity of this 
insight and try to understand the radically different roles played by Mg 2+ and Zn 2+ ions in 
DNA kink transition. The spontaneous curvature! of the DNA material used is anticipated 
to play a significant role and we model the DNA chain as a thin Kirchhoff elastic filament! 
with spontaneous curvatures, i.e., adjacent cross-sections of the filament are tilted a certain 
angle with respect to each other (in a mean-field sense). The filament will curve to a planar 
ring of curvature kq in its natural undistorted state, as what is observed in experiments.! 
Our study shows that a nonzero spontaneous curvature is required for the onset of the kink 
transitions observed in short DNA rings, in consistency with the insight of Han et al. 

We define a moving orthonormal coordinate system {d^s, t), d 2 (s, t), d 3 (s, t)} along the 
filament axis line, with d 3 being the tangential unit vector at arc length s. In general the 
filament configuration will change with time t,!0 so the direction vectors are also functions 
of t. The choice of the other two unit vectors di and d2 is arbitrary as long as di x d2 = d 3 ; 
in our case we define d± to point to the tilt direction of the filament cross-section at each 
arc length point s, for convenience. Then the free energy functional of such a filament is 
equivalent to the following form 

H = T I " Ko)2 + n " + r( ^ 3 " LJo)2 ) ds ((U) 

where E is the Young's modulus of the filament and / is the momentum of inertia of the 
filament cross-section (assumed to be circular) Jill T is a dimensionless parameter ranging 
between 2/3 and 1.II1JO f2 = yjjfijdj is called the twist vector, defined by the following 
equation 

d^ = nxd 4 , (2 = 1,2,3) 
here and after ( )' means s-derivative. We have chosen such a continuous and homogeneous 
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model partly because of the experimental fact that kinking locations are not (or only very 
weakly) correlated with base pair types.! 

In proposing Eq. (1) we have also take into account via the spontaneous twist rate uq 
the fact that ordinary DNA can also have a linking number deficiency of 5%. Ill For DNA 
circles with 168 bps (16 turns), it may well have 0.8 turns of initial twist in the DNA circle; 
for those with 126 bps (12 turns), it may well have 0.6 turns of initial twist. However, 
this (possible) spontaneous twisting tendency cannot explain the above mentioned kink 
deformations of short DNA. The reason is the following: For an elastic filament formed bv 
two chains interwinding around each other, such as DNA and actin, a well known modelli3li3 
concerning the effect of the fixed linking number shows that when the total twisting number 
variance ATw exceeds some threshold value Tw c = \/n 2 — 1A/C (n > 2), a planar circle 
will buckle to the n-th deformation mode (this relation is also rederived in Appendix B), 
where A and C are respectively the bending and twisting persistence length (for DNA A ~ 50 
nm and C ~ 75 nm.) If the kink deformations are indeed induced by this topological effect, 
then for the typical square (n = 4) polygonal kinked shape observed in the references^! to 
appear the total twisting number of the 168-bp DNA ring are required to deviate from its 
equilibrium value up to Tw c = 3 turns-I^HEl This value is almost four times the actual value 
of 0.8 turns. Thus it seems that the spontaneous twist will only play a neglectable role in 
the kink deformations of short DNA, and we feel it might be appropriate to focus on the 
macroscopic bending tendency of the filament by setting uo = in Eq. 1. The results derived 
based on this approximation also confirm this to be reasonable, as can be seen later. To be 
complete, we have also listed out the general result of model (1) in the case of u ^ atthe 
end of this paper (Appendix B), it reduces to the well known result in the referencest3ii3 
for the limiting case of DNA with intrinsic twist but no intrinsic bend (k q = 0). 

Recently Goriely, Tabor0 (and others) have suggested a powerful way to investigate 
on the dynamical properties of Kirchhoff thin filaments. We will use their procedure in 
studying the stabilities of model (1). Goriely and Tabor have worked out the stationary 
shape equations (SSEs) and the dynamical variational equations (DVEs)0 for the simplest 
case of a filament without any spontaneous curvature. Here, for our purpose, we first give 
the SSEs and DVEs for a general Kirchhoff filament. (Since referenced has demonstrated 
a way to get the DVEs we will not waste space in writing down the detailed calculations, 
but only notice here that there are some typographical errors in this reference.) After the 
general SSEs and DVEs are obtained, we then turn back to the special case of Eq. (1). 

The configurational free energy functional for a general Kirchhoff filament isETEJ 

ET r 

H = — J ((O, - K{f + (0 2 - Klf + T(Q 3 - Kl) 2 )ds, 

and the internal torque M is related to the deformation via the following constitutive equa- 
tion 

M = EI (Sit - Kf)d! + EI (Sk - iff)d 2 + EIT(0 3 - iff)d 3 , 

here (i = 1,2,3) is the spontaneous curvature along the dj direction. There are also 
internal forces F along the filament. We can perform a standard scaling as listed in Eq. 
(19) of reference^ to transform all the concerned quantities such as M and F into dimen- 
sionless forms. After this operation, we perform a basis perturbation operation suggested 
by Goriely and Tabor0 to the stationary configuration dj ^ and get that 
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d^df + eaxdf, (0.2) 
F = f^ + e(f^+f^axdi\ (0.3) 



M = (fil 0) - K?)dP + e (ajdf - e^of n<°> ), (0.4) 

where ( )(°) means value corresponding to the stationary state and ( p 1 ' its first-order 
correction, e is a small quantity. Insert Eqs. (2-4) into the Kirchhoff equations developed in 
references™^ and after a lengthy but straightforward calculation we can obtain the SSEs to 
be: 

(F<°))" = (A (0) dP + /f df + /f df )" = 0, (0.5) 

(fii 0) - k{)' - (of - i^of + r(of - = / 2 (0) , (o.6) 
(n (o) _ ^ 2 e y - r(n< 0) - ^)oS 0) + (of - /^)^ 0) = -/! {0) , (0.7) 

r(fif - KD' + i^fif - K^f = 0. (0.8) 

Eqs. (5-8) determine the stationary configurations and their corresponding internal force 
and torque (fif — Kf) distributions for a Kirchhoff filament with spontaneous curvatures. 
The stability of a stationary configuration is governed by the DVEs, which are the lengthy 
equations listed in Appendix A, Eqs. (24-29). 

The SSEs and DVEs derived above can be applied in studying the dynamical properties 
of any kind of Kirchhoff filaments, however in this paper we will only study a very special 
case, the elastic energy Eq. (1), with K\ = k and K 2 = K\ = 0. In this case, the SSEs 
Eqs. (5-8) demonstrate that the planar ring configuration is a stationary solution, with 

= k = const, flf = Qf ] = 0, 

A (0) = / 2 (0) = fP = o, (0.9) 

where k is the (actual) curvature of the ring. We now investigate on the dynamical stability 
of this ring shape. With Eq. (9), the DVEs Eqs. (24-29) reduce to 

(/i (1) )" = « 2 , (0.10) 

(/ 2 (1) r-2 K (/ 3 (1) y-^ 1) = -«i, (o.ii) 

(/i 1) ) /, -^/i 1) +2^(/ 2 (1) ) / = 0, (0.12) 

«i -/ 2 (1) =«i, (0-13) 
Tag + Yna' 2 + k§ol 2 — kkqol-?, = 2«3, (0-14) 
a 2 — TKa' 3 + (1 — T)K 2 a 2 — K a' 3 — K na 2 + /i = a 2 . (0.15) 

Redefining the perturbation parameters as (3\ = f[ , /3 2 = a 2 , /3 3 = tt3,/3 4 = f 2 \l3$ = 
fs,[3e = cui, the periodic solutions of this perturbation system Eqs. (10-15) are of the 
following form0 

Pj = e at (xje inKS + x*e~ inKS ) (j = 1, • • • , 6). 
Inserting this into Eqs. (10-15), we get, in matrix form, 

Lx = 0, (0.16) 
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where x = (xi, 22, £3, £4, x$, xq) t and L is defined as 

with being the 3x3 zero-matrix and 

/ -n 2 K 2 -a 2 \ 

Li = in(TK + K )K -Tn 2 K 2 - K K-2a 2 (0.17) 

\1 (—n 2 + l — T)K 2 — K K — a 2 —in(TK + K )K/ 

I -I -n 2 K 2 -a 2 \ 

L 2 = 2inn 2 -(u 2 + 1)k 2 (0.18) 



\ -(n 2 + 1)k 2 -2inn 2 



a 
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Eq. (16) has non-zero solution x if and only if the determinant of the matrix L equals to 

A L = A Ll • A L2 = 0, (0.19) 

where 

A Ll = -n 2 (n 2 - l)K 4 («g - (1 - T)k k - Tn 2 K 2 ) + a 2 ((2 + T)n 4 K 4 

-2(1 - I>V + 3n 2 K n 3 + Yn 2 K 2 + k q k + a 2 (2 + 4nV)), (0.20) 
A L2 = n 2 (n 2 - 1) V + kVV + 1) + kVV - l) 2 . (0.21) 

A^ 2 is related the excitation of {f^\ , «i} as can be inferred from Eq. (16) and the 
characteristics of the matrix L, and its value is not related to Kq, A^ 2 > for any n > 2. 
Then to satisfy Eq. (19) we need only to consider A^ which is related to the excitation of 
{/1 > a 2j a 3}- (We note that ^lAti = 1) = A i2 (n = 1) = only for a = 0, therefore the 
n = 1 modes are just soft modes" not important for the linear instability.) 

Solution of Al^ct) = with real positive a identify the unstable modes of {fi X \ 0.2, 03} 
(these modes will be unstable because the amplitudes of the small fluctuations will grow 
exponentially with time.) From Eq. (20) it is evident that this condition is equivalent to 
require 

g = k 2 - (1 - T)k k - Tn 2 K 2 > 0. (0.22) 

The general behavior of g is shown in Fig. 1. It clearly demonstrates that only when the 
actual curvature k of a ring is smaller than its spontaneous curvature k will it be possible 
for the ring to deform. This prediction is in agreement with the experiment of Han et a/.,Bi 
and the threshold curvature for the n-th mode to buckle is 



-i + r + A /(i-r) 2 + 4rn 2 

«c(n) = V -^— 2 Ko (n > 2). (0.23) 

For Kq — 0, g is always negative and no buckling process will occur. Therefore we can 
conclude that (i) k 7^ and (ii) k < k,q are the necessary conditions for the instability 
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of a ring (Fig. 1). Thus, our present work can qualitatively explain the novel phenomenon 
of DNA described in the introduction part: 126-bp DNA will not kink because no matter 
what ionic conditions its k is always > k c ; on the other hand, at some appropriate ionic 
conditions the curvature of a 168-bp DNA can become < k c , and kink deformation will be 
triggered.!^ When taking into account the possibility of nonzero ujq i n Eq. (1), a threshold 
condition similar to Eq. (23) is derived in Appendix B. 

At the initial stage of buckling, only {f[ , 0:2,0:3} will be excited, therefore the DNA 
configuration will deform out of the ring plane rather than deform in the same plane as the 
ring lies in. Thus in the present we can not tell whether the buckled shapes will actually 
evolve to the kinked ones observed in experiment .01 To know this, nonlinear analysis beyond 
the buckling point is needed and it is very much involved to perform. However, we believe 
that the linear analysis we employed here has correctly described the behaviors of DNA at 
the onset of kinking. 

Experimentsaia shows that Zn 2+ ions can enhance the instability of DNA rings, while 
Mg 2+ ions gives a negative effect. Comparing this with the theoretical analysis, it is reason- 
able for us to suggest that Zn 2+ and Mg 2+ ions will respectively increase and decrease the 
spontaneous curvature of the DNA chain. We are informed^! that that Zn 2+ mainly binds 
to DNA base pairs, and Mg 2+ interacts with the back-bone phosphate ions. It might be 
possible that the intercalation of Zn 2+ takes place mainly at the exposed side of the DNA 
ring hence causing an increase in its spontaneous curvature. The following analysis shows 
that the effect of Zn 2+ to the spontaneous curvature is very significant. 

The largest solution of Al 1 (ct) = changes with n at four fixed k/kq ratios are shown in 
Fig. 2. The possible deformation modes are obviously determined from these curves as they 
correspond to real positive values of a. The value of n corresponding to the peaks of these 
curves (n c ) are just the most observable modes Jll because these modes grow the fastest. At 
ImM ZnBr 2 it is shown experimentally that n c ~ 4,i therefore we estimate from the solid 
line of Fig. 2 that k/kq ~ 0.25. The value of k,q in Zn 2+ -free solutions is till not precisely 
known, but it must be higher than 27r/126 bp -1 .!'! On the other hand, in this condition 
168-bp is not kinked, so we must have k > k c (2), or in other words, k < 27r/90 bp -1 if we 
set T = 2/3 in Eq. (23). Therefore, when there is no Zn 2+ ions in the solution, 27r/126 bp -1 
< k < 27r/90 bp -1 , i.e., 0.55 < k/ko < 0.75, (the threshold case k/k = 0.55 is shown in 
the dot-dashed curve of Fig. 2.) Compare this estimate with the value 0.25 corresponding 
to ImM ZnBr 2 we can conclude that the addition of ImM Zn 2+ ions makes the value of the 
spontaneous curvature of DNA increase at least one times. 

This high efficiency of Zn 2+ make us believe it to be possible that Zn 2+ ions will induce 
spontaneous curvatures to DNAs even when they are originally linear and cause them to 
kink if they are loaded with stresses. This may be of vital biological significance. For 
example, it is known that transcription and replication of DNA occur with the participation 
of specific enzymes which contain the Zn 2+ ion at the active site. And it is also found that 
transcription proceeds simultaneously with conformational changes of the DNA chain.0 
More investigations on this respect are deserved. 

One of the author (Z. H.) appreciates valuable discussions with Dr. Yan Jie, Dr. Liu 
Quanhui and Dr. Zhao Wei. 
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Appendix A: The dynamical variation equations (DVEs) 

The dynamical stability of a certain stationary configuration of the general Kirchhoff 
filament are determined by the dynamical variation equations listed below in Eqs. (24- 
where ( ) means second order t-derivative: 

a 2 = [(a' 2 - a 3 nf ] + arf)' - {l?\a' 3 - a^P + 

- [(4 - arf + a 2 nf ] y - n 2 \a[ - a 2 n 3 0) + a 3 n {0) 
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+^\a' 2 -a 3 n? + ai n^)]f^ 



- 2(a' 3 - + ' - fPn? + /rni 0) ) 

+ A (1) " + 2/«'fi< 0) - 2/ 2 (1) '^ 0) - fi l) ((n { 3 y + (fif) 2 ) 

+ / 2 (1) (-^ 0) ' + ^°^S 0) ) + f?W + ^°^S 0) ), (0.24) 



-a l = [(q/ 3 - ai nr + a 2 nry - ft 2 Uj (4 - a 2 nr + a 3 n [ 2 U) ) 
+nf\a' 2 -a 3 n^ + ai nr)]A 0) 

- [(4 - a 2 Oj 0) + a^)' - nf\a' 2 - + arf) 

+ 2(4 - erf + a 2 nf>)(A< >' - /f Of + /f^ 0) ) 

- 2(4 - arf + a 3 ^ 0) )(/3 (0) ' - A ( °W + / 2 (0) fiS 0) ) 

+ A (1) " - + 2/«'n?> - / 2 (1) ((^! 0) ) 2 + (fif) 2 ) 

+ A (1) (-< + n£°>n<°>) + /W(nf + nW), (0.25) 

= [(4 - erf + erf)' - Q 3 °\a' 2 - + erf) 

rf(4-rf + erf)]/f 

- [{a! 2 - a 3 nf ] + arf)' - fi! 0) (4 - «i^ 0) + erf } ) 

rf(4-rf + rf)]A (0) 

+ 2(4 - erf + erf)(/ 2 (0) ' - /irf + A ( rf) 

- 2(a 2 - a 3 nf + erfXjf' - A ( rf + /|rf) 

+ /«" - 2A (1) '^ 0) + 2A (1) '^S 0) - / 3 (1) rf) 2 + (fif) 2 ) 

+ A (1) (rf' + ^rf) + / 2 (1 W + n?>n?>), (0.26) 



& x = (4 - erf + erf)' + (4 - a 3 n^ } + a^fxrrf - Ki) - «w) 

+ iT 2 e (4 - a x Slf + erf } ) + a 2 (r(^ 0) - K e 3 )' + Kl^f 1 - K e 2 ttf ] ) 

- a 3 rf - K°y - irf - K e 3 )n^ + (^S 0) - *rrf) - a (1) - « 3 A (0) , (0.27) 



7 



a 2 = (a' 2 - a 3 ^ 0) + «i^ 0) )' + (« 3 - «i^ 0) + a 2 ^ 0) )((nS 0) - K{) - TSlf) 

+ (ni 0) - r(^ 0) - K e 3 ))(a[ - a 2 ^ 0) + « 3 ^ 0) ) + « 3 ((^S 0) - K{)' - (fi<°> - K e 2 )nf ] 
+ T(^ 0) - ff|)n<°>)) - a 1 (r(fif - Kl)> + K{rt? - K^) + /« - a 3 / 2 (0) , (0.28) 

2a 3 = r(o / 3 - «i^ 0) + arf } )' - K e 2 (a[ - a 2 nf + a 3 fif) 

+ K{(a' 2 - a 3 fiS 0) + a&f ) + «i((^ 0) - K e 2 )' - Y(Q§ ] - + (Vl^ - K^ng*) 

- a 2 ((^ 0) - i^i)' - (fl< 0) - K e 2 )Qf ] + r(^ 0) - K^) + ai fl 0) + a 2 $\ (0.29) 
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Appendix B: The linear instability result for Eq. (1) with nonzero 

In the main text, we have focused our attention to the uq = case of medel Eq. (1). For 
the general case of nonzero uq, following the same procedure as discussed in the main text 
we can obtain the threshold condition for the n-th mode to become instable. It reads 

K 9 (n 2 — l)n 4 (n + l) 2 x 
(k 2 k (T - 1)(1 - n 2 ) + ««g(l - n 2 ) - T 2 u 2 (k + Tkti 2 ) + TK 3 n 2 (n 2 - 1)) = 0. (0.30) 

In the limiting case of ujq = (the filament has intrinsic bend but no intrinsic twist), Eq. 
(30) reduces to Eq. (23) of the main text. In the limiting case of /to = (then DNA has 
intrinsic twist but no intrinsic bend), Eq. (30) reduces to 

- lV nV - l) 2 (r 2 o; 2 - K 2 (n 2 - 1)) . (0.31) 

Eq. (31) is just the same result obtained in previous references for a chain with intrinsic 
twist jlHi! showing the correctness of this method. 
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FIGURES 



FIG. 1. The behavior of g [Eq. (22)] for n = 2 and V = 2/3. g is negative for k /k < 1 (the 
dotted part of the line), it becomes positive only when kq/k exceeds some threshold value higher 
than 1 (the solid part). g(tto/n = 1) = — T(n 2 — 1) and is negative for the important case of 
n > 2. Thus it is obvious that buckling process can take place only for k < kq. This prediction is 
confirmed by experiment. 1 ' 2 

FIG. 2. The relation between the largest solution of A^ 1 (cr) = [Eq. (20)] and n, at 
k/kq = 0.25 (the solid line), 0.30 (the dotted line), 0.40 (the dashed line), and 0.55 (the dot-dashed 
line) for the 168-bp DNA ring. We set k = vr/(0.34 * 168) for the 168-bp DNA 10 ' 12 ~ 15 and choose 
r = 2/3. 10 ' 15 The most dominant deformation mode determined by the peaks of these curves in- 
creases with the decreasing of k/kq, and the number of observable modes also increases as k/kq 
decreases. At k/kq = 0.25, the dominant mode is n = 4 and the modes n = 2, 3, 4 can be observed; 
while at k/kq = 0.55, no deformation mode with n > 2 will be excited. This prediction is in very 
close agreement with experimental observations. 2 



11 



